Homotopy Theory of Topological Defects in Spinor Condensates 
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We investigate the topological defects in atomic spin-1 and spin-2 Bose-Einstein condensates 
by applying the homotopy group theory. With this rigorous approach we clarify the previously 
controversial identification of symmetry groups and order parameter spaces for the spin-1 case, and 
show that the spin-2 case provides a rare example of a physical system with non-Abelian line defects, 
and the possibility to have winding numbers of 1/3 and its multiples. 
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The all-optical trapping of Bose-Einstein condensates 
(BEC) [l| has opened up a new direction in the study 
of dilute atomic gases, i.e., the spinor condensates with 
degenerate internal degrees of freedom of the hypcrfinc 
spin F. For alkali atoms with F = 1, both experiments 
and theories have shown two possible kinds of spin cor- 
relations in the atom species, namely ferromagnetic (e.g. 
87 Rb) or antiferromagnetic (e.g. 23 Na). With the exper- 
imental success of condensing alkali bosons with F > 1 
such as 85 Rb and 133 Cs 2], and the unusual stability 
of the F — 2 state (against spin-exchange) in 87 Rb [3|, 
one expects that defects with much richer structure can 
be created in the future. A remarkable feature here is 
that not only the gauge symmetry U(l) but also the spin 
symmetry SO (3) are involved, a situation similar to su- 
perfiuid 3 He where three different continuous symmetries 
(orbital, spin and gauge) are broken either independently 
or in a connected fashion Q. 

Topological defects and excitations in the spinor BECs 
have been studied theoretically by several groups 0, 0, 
0, H E3, EH Most of the work on this subject is 
based on the original identifications of the order param- 
eter spaces by Ho After the original studies it was 
also claimed by Zhou that a discrete symmetry of Z-j 
type was missed in the case of antiferromagnetic spin- 
1 condensate Lltj and therefore the topological defects 
would manifest totally different structures. In this Letter 
we present a rigorous topological study that both solves 
this spin-1 controversy, and reveals interesting aspects of 
spin-2 systems. The phases of spin-2 spinor condensates 
are characterized by a pair of order parameters (F) and 
|0| describing the ferromagnetic order and the formation 
of singlet pairs, respectively 0, 0, Q|. It turns out 
that for the so called cyclic phase the fundamental group 
that determines the nature of possible stable topological 
defects is non-Abelian. The only known physical exam- 
ple of such a system so far has been the biaxial nematic 
liquid crystal. 

Outline of the homotopy theory of defects: We sketch 
out the procedure which has been widely used in the 
study of topological defects in ordered media such as 



liquid crystals, superfluid 3 He and heavy-fermion super- 
conductors [15j . The central feature of the classification 
scheme of the defects emerges from examining the map- 
pings of closed curves in physical space into the order- 
parameter space (OPS). The order parameter of a spinor 
BEC, 'F(r) = */ n(r)£(r), where n(r) is the density and 
£(r) is a normalized spinor £t ( r ) . £( r ) = 1, has associated 
with it a group of transformations G, i.e., all spinors re- 
lated to each other by gauge transformation e and spin 
rotations e~ tF * a e~' lFyl3 e~' tF *~ l are energetically degener- 
ate in zero external magnetic field, where {a, (3, 7) are 
the Euler angles. The set of all transformations in G 
that leave the reference order parameter (chosen arbi- 
trarily but thereafter fixed) unchanged is known as the 
isotropy group H . The OPS can then be taken to be the 
space of cosets of H in G: M = G / 'H . In terms of broken 
symmetry, the fact that the ordering breaks the under- 
lying symmetry is expressed in the fact that H is only a 
subgroup of G. Let G be a connected, simply connected 
continuous group and Ho be the set of points in H that 
can be connected to the identity by a continuous path 
lying entirely in H. The fundamental theorems assert 
that 



7n(M) = H/H , ir 2 (M) = m(H ) 



(1) 



where the fundamental (second) homotopy group 7ri(7r2) 
characterizes the singular line (point) defects in three- 
dimensional space. 

For the spinor condensate it seems natural to iden- 
tify the underlying symmetry group as U(l) x SO(3), the 
groups in the direct product representing the gauge and 
spin degrees of freedom respectively. This group is not 
simply connected, i.e., 7n(U(l) x SO(3)) 7^ 0. To apply 
the theorems, however, it is essential that one chooses 
the group G to be simply connected (i.e., -K\{G) = 0). 
This can always be done because any continuous group 
can be imbedded in its universal covering group. We 
proceed by specifying the symmetry group as R x SU(2), 
with the group of real numbers R representing any trans- 
lation £ (-co, +00) in the phase of the condensate. 
For F = 1, we use the 3D representation of the group 
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SU(2) in order to obtain the isotropy group, e.g., a 
rotation u(z, a) around axis z by angle a takes the 
form of a diagonal matrix Diag (e~ M , l,e m ). The ele- 
ment — u(z, a) is represented by the same matrix, i.e., 
u(z, a + 2ir) — —u(z, a) while u(z, a + 4ir) — u(z, a). 

Calculation of the homotopy groups: There are two 
possible ground states in F = 1 case. For the ferromag- 
netic state, the isotropy group H is constructed by the set 
of transformations which leave the reference order param- 
eter (1,0, 0) T invariant. From the ground state spinor 

£ = e i(0-TO (e~ ia cos 2 f , V2 cos | sin | , e ia sin 2 §] we 
know immediately that the angles should satisfy 



/3 = O,0-or-7 = 2nn 



(2) 



with n an integer. The elements in group H are 
the combination of a translational part and a ro- 
tational part H = {{9,u(z,9)),(9,u(z,9 + 27r))} = 
{(9, ±it(z, &))}■ Evidently this group includes two dis- 
connected components — the connected component of the 
identity Hq = {(9,u(z,9))} is isomorphic to R. The 
group H/Hq is isomorphic to the integers modulo 2, i.e., 
Z 2 . The second homotopy group tt 2 is trivial and we 
arrive at the same result as that in Ref. |9( 

tti(M) = Z 2 , tt 2 (M) = 0, (spin-1 FM state). (3) 

A ferromagnetic spin-1 condensate may have therefore 
only singular vortices with winding number one while 
the point-like defects are topologically unstable. 

The polar state emerges if the atoms in the conden- 
sate interact anti-ferromagnetically. In the ground state 
C = e i6 /V2 (-e~ ia sin/3, v^cos/?, e M sin/?) T , the refer- 
ence parameter (0,1, 0) T is left invariant for just those 
elements with 

(3 = 0,9 = 2mr or (3 = tt, 9 = (2n + 1) tt. (4) 

Thus the isotropy group H must now be expanded to 
include not only transformations in which both the rota- 
tion and the translation leave the spinor unchanged, but 
also those in which the rotation takes the reference spinor 
(0, 1, 0) T to (0, —1, 0) T and the translation takes it back, 
i.e., a tt rotation about arbitrary axis perpendicular to z 
combined with a tt translation in 9 (or any odd multiples 
of tt). The full isotropy group is the union of these two 
sets, H — {(2rnr, u(z, a)) , ((2n + 1)tt, gu(z, a))} where 
g = u(y,7r). There are infinitely many discrete compo- 
nents in H, while the connected component of the iden- 
tity H = {(0, u(z, a))} is isomorphic to U(l). The ele- 
ments with an even translational parity are of the form 
(2mr, I)Hq, and those with an odd parity are of the form 
((2n + 1)tt, g) H . The group H/H is therefore isomor- 
phic to the group of integers Z through the isomorphism 
((2n + j)tt, gj) H i— > 2n + j for j — 0,1. We recover 
the conclusion that line and point defects in spin-1 polar 



state can be classified by integer winding numbers, 

m(M) = Z, tt 2 (M) = Z, (spin-1 Polar state). (5) 

Thus the Z 2 term does not appear in the homotopy 
grou p. A s we shall see the identification of the OPS in 
Ref. [l0( is also incorrect. 

Spin-2 Bose condensate: We next apply the same ap- 
proach to the BEC of spin-2 bosons. There are three 
possible phases for a spin-2 spinor Bose condensate, i.e., 
one more compared to the spin-1 case. This extra phase 
comes from the additional interaction parameter describ- 
ing the singlet formation, and is referred to as the cyclic 
phase [l2l H?3 . Il^j . The defects which may be cre- 
ated in spin-2 condensate exhibit even more elaborate 
structures due to quantum correlations among bosons. 
For F = 2 we have to use the 5D representation of 
SU(2), e.g., the rotation u{z,a) is represented by ma- 
trix Diag (e~ 2ia , e~ la , 1, e iQ , e 2la ). The calculations of 
the degenerate family of the ground state spinors and 
the corresponding homotopy groups are straightforward 
and details will be presented elsewhere j!6| . Here we only 
pick up some interesting features in our results, focusing 
on the symmetry properties of the defects in comparison 
with those in other ordered media. 

Equating the general expression for the ground state 
spinor of the ferromagnetic state F 



-1(0-27) 



/ e" 2M cos 4 f \ 
e~ la sin (3 cos 2 § 

^sin 2 /3 
e ia sin (3 sin 2 f 
\ e 2m sin 4 § J 



(6) 



with the reference spinor (1,0,0, 0,0) T leads to the re- 
quirement for the isotropy group H 



(3 = 0, 9 - 2a - 2 7 = 2mr. 



(7) 



We see that taking n — 0, 1,2,3 is enough for all pos- 
sible transformations, with the translational part being 
arbitrary and the rotational part containing the rota- 
tions around axis z by 9/2, 9/2 + tt, 9/2 + 2vr, 9/2 + 3tt 
respectively. Hence the group H is composed of 4 pieces 
H = {(9, u(z, mr + 9/2))}. Here it is important to show 
that the 4 components are not connected: there does not 
exist a continuous path in H which connects one compo- 
nent to another, though the rotational parts themselves 
are connected. The connected component of the identity 
Ha = {(9, u(z, 9/2))} is again isomorphic to R. If we de- 
fine an element g of the group R x SU(2) by (0, u(z, tt)), 
we see that the quotient group H/Hq has the same struc- 
ture as the cyclic group of order 4, i.e., {e, g, g 2 , g 3 } and 
we conclude that 



771 (M) 



7r 2 (M) = 0, (spin-2 F state). (8) 



It is interesting to check how the group Z<± character- 
izes vortices for state F. In spin-1 case there is only 
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one topologically stable line defect, that is, a vortex with 
winding number one. Equation JSJ shows that there 
are 3 stable vortices for spin-2 condensates. We can 
set 6 — 27 = 2mip, —a — m<p, (3 — 7rf in the ground 
state for F state, Eq. ©, which leads to a family of 
spinor states parametrized by a parameter t between 
and 1. Here m > is an integer, tp is the azimuthal 
angle. When t evolves from to 1, the 4nup vortex 
state C(t = 0) = (e i4mv ,0,0,0,0) T evolves continu- 
ously to the vortex free state £(£ = 1) = (0, 0, 0, 0, 1) T . 
This shows that vortices with winding number Am are 
topologically unstable. Similarly, by multiplying factors 
e tkip (k = 1,2,3) one obtains the following correspon- 
dences 

e i(4m+k) v ^ Qj Qj jT ^ eikv ^ Qj Q) j_jT ^ 

i.e., the vortices with winding numbers Am+k may evolve 
into vortices with winding numbers k, respectively. There 
are thus 3 classes of topologically stable line defects. To- 
gether with the uniform state, they form the fundamental 
group Z4. Straightforwardly for condensates with spin 
F, the fundamental group 7i"i(M) = Z 2 f characterizes 
(2F — 1) classes of stable line defects. 

Non-Abelian homotopy groups: Media with non- 
Abelian fundamental groups are especially interesting 
from the topological point of view. The only illustra- 
tive example in ordered media so far have been biaxial 
nematic liquid crystals 01- ^ ^ i s taken to be SO(3) 
then the isotropy group H is the four-element group, D 2 , 
the symmetry group of a rectangular box. If, however, 
we take G to be SU(2), then H is expanded to the non- 
Abelian quaternion group Q (known as the lift or double 
group) with eight elements, the multiplication table of 
which has been verified experimentally |l8j . 

We have found that the cyclic state C provides an- 
other physically realistic example in which the funda- 
mental group is non-commutative. The reference spinor 

\ (e 1 *, 0, v2j 0, — e~ 1 ^) is left invariant by the elements 
in the isotropy group 

H = {±I,±a,±b,±c, (10) 
±d, ±e,±f,±g, 
±d 2 ,±e 2 1 ±f 2 ,±g 2 }. 

The spin rotations a — u(z, 7r), b = u(y, tt)u(z, <f> + tt/2) 
and c = ba satisfy a 2 — b 2 — c 2 — —I, while 
d = u(z, 7t/4 + <j)/2)u(y, n/2)u(z, n/4 - cj>/2), e = -da, 
f = —ad and g — —ada satisfy d 3 — e 3 = f 3 = g 3 = —I. 
Each element in the first, second, third row is associated 
with an additional phase change 2nn, 2ir/3 + 2mr, 47r/3 + 
2n7r respectively. It is a discrete group, and H consists 
of the identity (0, 1) only. The fundamental theorems 
identify that 

7Ti (M) = H, tt 2 (M) = 0, (spin-2 C state). (11) 



The elements in the fundamental group are non- 
commutative, for example ab = —c =/= ba. The criterion 
for the topological equivalence of defects applies in the 
most general case in terms of conjugacy classes of the 
fundamental group. It is thus necessary to classify the 
group into the following conjugacy classes: {/}„, {—/}„, 
{±a,±6,±c}„, {d,e,f,g} n+1/3 , {-d,-e,-f,-g} n+1/3 , 
{d 2 ,e 2 , f 2 , 3 2 }„+2/3, {-d 2 , -e 2 , -f 2 , -g 2 } n +2/3 with the 
subscripts standing for the winding numbers of the de- 
fects. Physically this indicates the feasibility of creating 
not only vortices with any integer winding number but 
also fractional quantum vortices. The homology theory 
fli^ l assembles the conjugacy classes further into 3 sets for 
each n, in which the defects are labeled by the winding 
numbers n,n + 1/3, n + 2/3 respectively |l6(. 

Interesting features of this non-Abelian fundamental 
group include the topological instability of the defects 
and their interaction, i.e., entanglement when two of 
them are brought to cross with each other [l^. Two de- 
fects in the same conjugacy classes can be continuously 
converted into one another by local surgery. In addition, 
the defects in the same homology set can be transformed 
into each other via a catalyzation process consisting of 
splitting a line singularity into two and recombining them 
beyond a third one 0] . 

Like the quaternion group Q for biaxial nematics, the 
fundamental group (|ll)f) is the lift of a point group in 
R x SU(2). To find the remaining discrete symmetry 
group for the cyclic state, and, in addition, to clarify the 
controversial identification of the OPS for spin-1 case, 
in the remaining of this paper we turn to describe the 
system in terms of rotations in SO (3), e.g., two elements 
±u(z, a) in SU(2) are mapped into one R(z, a) in SO(3) 
with R(z, a + 2ir) = R(z, a). 

Order Parameter Spaces: The OPS for F = 1 polar 
state was identified as U(l) x S 2 in Refs. 

An ex- 
tra Z2 symmetry was claimed in Ref. fu\ so the author 
concluded the OPS as U(l) x S 2 /Z 2 - Here we show that 
previous studies are incorrect. Taking the group G as 
U(1) G x SO(3) s where the subscripts stand for the gauge 
and spin symmetries respectively, we see that the isotropy 
group H consists of two separate parts, { (e*°, R(z, a)) } 
and {(e l7r , R(y, n)R(z, a))}. The rotations in the first 
part constitute the group SO (2), while the elements in 
the second part are just those in the group 0(2) but 
not in SO(2) with determinants — 1. The combination 
of these two parts gives the full isotropy group as 0(2) 
where both gauge and spin symmetries are involved. The 
OPS is the quotient G/H = (U(1) G x SO(3) s ) /0(2) G+S 
and here it is not possible to apply the fundamental the- 
orem for G is not any more simply connected. 

For the ferromagnetic state the group H may be ob- 
tained if one notices that the 2tt difference in the rota- 
tional angle does not give another component as it did 
in the case of SU(2). We have H = {{e ie , R(z , 8))} 
which is isomorphic to U(1) G+S . This means that 
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there is a remaining symmetry U(l) in the symme- 
try broken system. The OPS is thus factorized as 
(U(1) G x SO(3) s ) /U(1) G+S = SO(3) s+G . 

The discrete symmetry group of defects in the spin- 
2 cyclic state C can be shown to be isomorphic to the 
tetrahedral group T. We continue to represent G as 
U(1) G x SO(3) s . The isotropy group ifTUI) is shrunk to a 
group of 12 elements if one understands the rotation in 
the sense of SO(3) (i.e., a — R(z,nj), 

H = {I, a, b, c, ed, ee, ef, eg, e 2 d 2 , e 2 e 2 , e 2 f , e 2 g 2 }, (12) 

where e = exp(2-7rz/3) comes from the gauge transfor- 
mation and ed, for instance, is an abbreviation for the 
element (e,d). Three 2-fold rotational axes are z and 
2 lines in xy plane perpendicular to each other (which 
lie on axes x and y if we choose the arbitrary phase 
4> = 7r/2). The elements ed,ee,ef,eg are four 3-fold 
axes. The symmetries remaining in the symmetry bro- 
ken states for biaxial nematics and spin-2 cyclic state are 
shown in Figure 1. The OPS for state C can be identified 
as (U(1) G x SO(3) s ) /T G+S . 
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FIG. 1: Symmetries of the defects in biaxial nematics (-D2) 
and cyclic state C in spin-2 condensate (T). The dot in the 
center of the rectangle stands for axis z. The dashed lines 
represent 2-fold axes, except that with a triangle for 3- fold 
axis. 

It should be noted that an applied magnetic held 
B changes the defect structure severely by reduc- 
ing the degenerate family of the spinor. We take 
again the cyclic state as an example. The sym- 
metry group in this case is U(l) x SO (2) because 
the magnetic held chooses its direction automat- 
ically as the quantization axis. From the spinor 

/ \ T 

\ f (1 + p)e^\0, ^2 - 2pVW ,1+ ^/ 2 ,0, (-1 +p)e i <M , 

where 0i,2 are two arbitrary phases and p ~ B, we 
easily see the possibility to create vortices in any of 
the three nonzero components with winding number for 
i— th component rij conhned by n\ + TI5 = 2713. 



In summary, we have determined the nature of the 
topological defects in spin-1 and spin-2 condensates. The 
order parameter spaces are identified as the spaces of the 
cosct of the isotropy group H in the transformation group 
G. Topologically stable vortices with winding numbers 
larger than unity may be created in the ferromagnetic 
state for condensates with F > 1, up to the value (2F — 
1). The line defects in the spin-2 cyclic state C exhibit 
non-commutative features, resulting e.g. in line defects 
with winding numbers of 1/3 and its multiples. 
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